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Abstract 

This article deals with the issues of global-in-time existence and asymptotic 
analysis of a fluid-particle interaction model in the so-called bubbling regime. 
The mixture occupies the physical space f2 C R'^ which may be unbounded. 
The system under investigation describes the evolution of particles dispersed in a 
viscous compressible fluid and is expressed by the conservation of fluid mass, the 
balance of momentum and the balance of particle density often referred as the 
Smoluchowski equation. The coupling between the dispersed and dense phases is 
obtained through the drag forces that the fluid and the particles exert mutually 
by action-reaction principle. We show that solutions exist globally in time under 
reasonable physical assumptions on the initial data, the physical domain, and 
the external potential. Furthermore, we prove the large-time stabilization of the 
system towards a unique stationary state fully determined by the masses of the 
initial density of particles and fluid and the external potential. 



1 Introduction 

Fluid-particle interaction arise in many practical applications in biotechnology medi- 
cine [4], reactive gases formed by fuel droplets in combustion [1, 2, 25, 26], recycling 
and mineral processes [5, 2.3], and atmospheric pollution [24]. Aerosols and sprays 
can be modelled by fluid-particle type interactions in which bubbles of suspended 
substances are seen as solid particles. Two-phase flow hydrodynamic models have 
also been proposed [(>]. 

Here, we focus on a particular system derived by formal asymptotics from a meso- 
scopic description. This is based on a kinetic equation for the particle distribution 
of Fokker-Planck type coupled to fluid equations. Different macroscopic equations 
can be obtained as scaling limits, see [7] for a complete description of the modelling 
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issues. In these models, the fluid is either incompressible [15, 16] or compressible 
[7] . The coupling between the kinetic and the fluid equations is obtained through the 
friction forces that the fluid and the particles exert mutually. The friction force is 
assumed to follow Stokes law and thus is proportional to the relative velocity vector, 
i.e., is proportional to the fluctuations of the microscopic velocity ^ € K'^ around the 
fluid velocity field u. More precisely, the cloud of particles is described by its distri- 
bution function /^(t,^,^) on phase space, which is the solution to the dimensionlcss 
Vlasov-Fokker-Planck equation 



9Je + ^(C-V,/,-V,$-V4/,) =- diV5((C-yi«e)/ + V5/e). (1.1) 



Here, e > is a dimensionlcss parameter and we have a a drag force independent 
of the fluid density g^, but proportional to the relative velocity of the fluid and the 
particles given by ^ — u^{t,x). The right hand-side of the moment equation in the 
Navicr- Stokes system takes into account the action of the cloud of particles on the 
fluid through the forcing term 



The density of the particles r}^{t,x) is related to the probability distribution function 
fe{t,x,^) through the relation 



The well-posedness of this kinetic-hyperbolic coupled fluid-particle system has 
been addressed in [19] in the case of compressible models for the fluid equations. There 
are two different scaling limits for this model, the so-called: bubbling and flowing 
regimes. They correspond to the diffusive approximation of the kinetic equation, 
bubbling regime, written in (1.1), and the strong drag force and strong Brownian 
motion for the flowing regime. This last regime has been studied in [2i i] where it 
is obtained rigorously as the limit from the mesoscopic description for local in time 
solutions and initial data bounded away from zero for the densities. We also refer to 
[8] for asymptotic preserving numerical schemes in relation to these scaling limits. 

In all the above mentioned studies, the viscosity of the fluid was neglected although 
it is the source of the drag forces. The viscosity is present in the dimensionless systems 
although negligible as noted in [~, Remark 3]. In this work, we will deal with the 
resulting formal macroscopic fluid-particle system obtained through the scaling limit 
in (1.1) as e by the standard Hilbcrt-cxpansion procedure. In this scaling limit, 
particles are supposed to have a negligible density with respect to the fluid, and thus, 
due to buoyancy effects, they will typically move upwards in a system under gravity, 
from where it bears the name of " bubbling" . This situation is typically complemented 
with no-flux boundary conditions in a bounded domain. More generally, we can ask 
ourselves under what conditions on the external potential in unbounded domains we 
can assert convergence towards stationary integrable states. We refer to [7, 8] for a 
detail account of the physical meaning and validity of the scaling limits. 

Summarizing, the state of such flows in this macroscopic description is, in general, 
characterized by the variables: the total mass density g{t, x), the velocity field u(t, a;), 
as well as the density of particles in the mixture rj{t,x), depending on the time t G 
(0, T) and the Eulerian spatial coordinate a; G 17 C M'^. 
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In this section, we present the primitive conservation equations governing fluid- 
particle flows in the bubbling regime. These equations express the conservation of 
mass, the balance of momentum, and the balance of particle densities often referred 
as the Smoluchowski equation: 

dtg + div,{gu)=0, (1.2) 



dt{gu) + div^(£.u (g}u) + V^(p(^?) + 7]) - fiAu - AV^^div^u = -(ry + /3g)W^^, (1.3) 



dtr] + div^(77(u - V^$)) - A77 = 0. 



(1.4) 



Here, p denotes the pressure p{g) — ag'' , a > 0, 7 > 1, /3 7^ 0, and (f> denotes the 
external potential (typically incorporating gravity and boyancy). 

In this paper, we require the potential to satisfy suitable confinement conditions 
(HC) (see Section 2), which docs not limit the physical relevance of our results. The 
viscosity parameters fi > and A + |/i > are nonnegative constants and /? > if 
is unbounded. 

Another macroscopic effect is that the total pressure function in the momentum 
equation depends on both the particle and the fluid densities p{g) + rj. We impose 
the no-slip boundary condition for the velocity vector leading to no-flux for the fluid 
density through the boundaries and the no-flux condition for the particle density 



u\on = ^xV ■ ^ + V^x^ -v^O on (0, T) x dn, 



(1.5) 



with ly denoting the outer normal vector to the boundary d^l. Our problem is supple- 
mented with the initial data {go, mo, 770} such that 



£>(0,x) = gae L'<{n)r]Ll{n), 
{gu){0,x) = mo e L^n)r\L\n), 

ry(0,a;) = r/o G L^{n)nL\{n). 



(1.6) 



Motivated by the stability arguments in [7], the numerical investigation presented in 
[8] , as well as a number of studies on numerical experiments and scale analysis on the 
proposed model (see [3]), we investigate the issues of global existence and asymptotic 
analysis for fluid-particle interaction flows providing a rigorous mathematical theory 
based on the principle of balance laws. The total energy of the mixture is given by 



Eirj, g, u)it) / \^g{t)\nit)\' + ^g\t) + (rj \ogrj)it) + {^g + 7?)(i)$ 

Jo L ^ 7-1 



dx (1.7) 



At the formal level, the total energy can be viewed as a Lyapunov function satisfying 
the energy inequality 



dE_ 
Itt 



dx < 0. 



fi\V^uf + A|div^u|2 + |2V:,v^ + 07V,.$|2 

Therefore, it is reasonable to anticipate that, at least for some sequences t„ — > cxd 
r]{tn)^ris, g(tn) ^ gs, £'u(t„) -> 0, 



(1.8) 
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where rjs , gs satisfy the stationary problem 

^xipiQs) + vs) = -ivs + /^Psjv^.* on n. 

The energy estimate written in the form 

Eit)+ f i\\^xu\\l. + \\dW.Mh)dt+ [ / |2V,y^+V^V,$p dxdt <i;(0), 
Jo Jo Jn 

now implies that 

|2V,V^+ V^V,$|2 =0. (1.9) 

The aim of this paper is to show that, in fact, any weak solution converges to a fixed 
stationary state as time goes to infinity, more precisely, 

g{t) — > Qs strongly in L''{n), 

esssup / £i(r)|u(T)p dx — > 0, 
T>t Jn 

•q{t) 7]s strongly in L^{rt), 

as t — >■ oo under the confinement hypothesis on the domain fl and the external 
potential $ given in (HC) (cf. Definition 2.1). 

Indeed, it can be shown that the sequences {Qm^n) of the time shifts defined as 

Qn{t,x) -.^ g{t + Tn,x), T„->00, 

r]n{t,x) :=r]{t + Tn,x), T„ CXJ, 
contain subsequences, denoted by the same index w.l.o.g., such that 

Qn -> Qs strongly in iL((0, 1) x fi), 

and 

Vn — > rjs strongly in ((0, T); L^^)(f2) for some pi,P2 > 1, 
where (ps,?7s) solve the stationary problem 

It is worth noting that the confinement hypothesis is both necessary and sufficient for 
the stationary problem (1.10) to admit a unique solution {gs,T]s) (Section 4). 
The main ingredients of our approach can be formulated as follows: 

• A suitable weak formulation of the underlying physical principles is introduced. 
Physically grounded hypotheses are imposed on the system: The mixture oc- 
cupies the physical space C M"^. The boundary conditions are chosen in such 
a way so that the dissipation of energy is guaranteed, whereas the pressure of 
the mixture takes into account both the density of the fluid and the density of 
particles. 

• A priori estimates are established based solely on the boundedness of the initial 
energy and entropy of the system. 
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• A suitable approximating scheme is introduced for the construction of the solu- 
tion based on a two level approximating procedure: The first level involves an 
artificial pressure approximation, whereas the second level approximation em- 
ploys a time discretization scheme. The sequence of approximate solutions is 
constructed with the aid of a fixed point argument coupling the time-discretized 
compressible isentropic Navier-Stokcs equations to a discretization in time of the 
equation for rj. 

• Physically grounded hypotheses are imposed on the domain D, and the external 
potential $ (confinement hypotheses (HC)). The analysis in the present article 
treats both the case of a bounded physical domain 57 as well as the case of 
an unbounded domain. We remark that the confinement hypothesis (HC) on 
(57, $) plays a crucial role in providing control of the negative contribution of 
the physical entropy t] log rj in the free-energy bounds for unbounded domains. 

• High integrability properties for the density need to be established for the limit 
passage in the family of approximate solutions and in particular in taking the 
vanishing artificial pressure limit. We remark that in the present context, the 
potential $ is not integrable on unbounded domains. To deal with this new 
difficulty we employ the Fourier multipliers in the spirit of [If, f8], while taking 
into consideration the new features of our problem. 

• We remark that both the total fluid mass and the total particle mass are con- 
stants of motion. In particular, we are able to conserve the total masses also in 
the large-time limit allowing us to uniquely determine the long time asymptotics 
(cf [14]). 

The paper is organized as follows. In Section 2 we collect all the necessary hypothe- 
ses imposed on the external potential (confinement hypotheses (HC)), we present the 
notion of free energy solutions and the main results of this article. Section 3 is devoted 
to the proof of the global existence of weak solutions (Theorem 2.1). First, a suitable 
approximation scheme based on an artificial pressure approximation and on a time 
discretization scheme is introduced. The remaining section is devoted to the limit 
passage in the family of approximate solutions. The most delicate part of the anal- 
ysis concerned with the vanishing artificial pressure limit is presented in Subsection 
3.5. The large time asymptotic analysis is described in Section 4. 

2 Free-energy solutions and main results 

In this work, we analyse the existence and large-time asymptotics of certain type 
of weak solutions to the two-phase flow problem (1.2)-(1.4) coupled with no-fiux 
boundary conditions (1.5) in two different geometrical constraints of interest in the 
applications: for bounded domains and for unbounded domains under confinement 
conditions due to the external potential. We will collect all assumptions concerning 
the geometry Q and the external potential $ under the generic name of confinement 
conditions. Let us remark that the external potential $ is always defined up to a 
constant. Therefore, for bounded from below external potentials we can always 
assume without loss of generality, by adding a suitable constant, that 

inf $(a;) = 0. (2.1) 
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Definition 2.1 Given a domain O S C'^''^ , v > Q, VI G M^, and given a hounded from 
below external potential $ : O — > satisfying (2.1), we will say that ($7, $) verifies 
the confinement hypotheses (HC) for the two-phase flow system (1.2)-(1.4) coupled 
with no- flux boundary conditions (1-5) whenever: 

(HC-Bounded) // is bounded, <& is bounded and Lipschitz continuous in 
and the sub-level sets [$ < k] are connected in $7 for any k > 0. 

(HC-Unbounded) // is unbounded, we assume that $ e w^^^in), I3>0, 
the sub-level sets [$ < k] are connected in Q for any k > 0, 

and 

|A$(a:;)|<ci|V,$(x)|<C2$(a;), > i?, (2.2) 
for some large R > 0. 



Remark 2.1 The confinement assumption (HC) has physical relevance in our set- 
ting as it is verified for several domains with $ being the gravitational potential. 
For instance, 

1. when Vl = {x <eM? \ {xi,X2) € [a,^]^, 2^3 G [0,^^]} and (^{x) = gx^, where 

^ Qp 

2. when Q, = {x {xi,X2) G [a, ^]^, a^s > 0} and <^{x) = gx^, where P = 
1 ~ ^ and qf < Qp- 

3. when = R'^ \ B{0,R) and $(a;) = g\x\, where B{0,R) is the ball centered at 
the origin with radius R and /3 > 0. 

Here, Qp and gp are the typical mass density of fluid and particles, respectively. 
Remark that 1. corresponds to the standard bubbling case (see [7j) in which particles 
move upwards due to buoyancy. 

Let us now specify the type of weak solutions for the two-phase flow system (1.2)- 
(1.4) we will be dealing with. 

Definition 2.2 Let us assume that (^7, $) satisfy the confinement hypotheses (HC), 
we say that {g,u^r]} is a free-energy solution of problem (1.2)-(1.4) supplemented 
with boundary data for which (1.5) holds and initial data {go,inQ,rjo} satisfying (1.6) 
provided that the following hold: 

• p > represents a renormalized solution of equation (1.2) on (0,cx)) x VL: For 
any test function ip G 2?([0, T) x 17), any T > 0, and any b such that 

6 e n C[0, 00), B{g) = B{1) + ^ dz, 

Ji z 

the following integral identity holds: 
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/ / (B{g)dtv^ + B{g)u ■ V,(/j - b{g)dw,u^) dx dt = - f B{goMQ, •) dx (2.3) 
Jo Jn ^ ' Jn 

• The balance of momentum holds in distributional sense, namely 

J J (^gu ■ dtv + gu ig) u : V^v + {p{g) + 77) div^jV^ d.T dt = (2.4) 

/ / (m^suV^jV + Adivj^udivj^v — {rj + Pg)\/x^ ■ v dx dt — / mo • v(0, •) dx, 
Jo Jn Jn 

for any test function v G I'([0, T); T>(^; R^)) and any T > satisfying <^|ao = 0. 

All quantities appearing in (2.4) are supposed to be at least integrable. In par- 
ticular, the velocity field u belongs to the space L^(0,T;W^'^{fl; R^)), therefore 
it is legitimate to require u to satisfy the boundary conditions (1-5) in the sense 
of traces. 

• rj >0 is a weak solution of (1.4). That is, the integral identity 

/ / r]dt(p + r]u-'Vj;(p-r]'^x^-'^xyy-'^xV'^xyy dxdt = - / r]o(p{0,-)dx (2.5) 
Jo Jn Jn 

is satisfied for test functions (ys e I?([0, T) x f2) and any T > 0. 

All quantities appear in (2.5) must be at least integrable on (0,T) x fi. In 

particular, 77 belongs L^{[Q,T];L^{Vt))CM^{0,T;W^^i{^)). 

• Given the total free- energy of the system by 

E{g,u,ri){t) Qg|u|2 + -^pT + 77log77+(^e + ?/)$^ dx, 

then E{g,\i,vi){t) is finite and bounded by the initial energy of the system, i.e., 
E{g,\i,ff){t) < E{gQ,UQ,ri(j) a.e. t > 0. Moreover, the following free energy- 
dissipation inequality holds 

[ / (M|V,u|2 + A|div,up + |2V,V^+V^V,$|2) d.Tdt<i;(f?o,uo,77o). (2.6) 
Jo Jn 

Now, we have all ingredients to state the main results of this work. 



Theorem 2.1 (Global Existence) Let us assume that (1^,$) satisfy the confine- 
ment hypotheses (HC). Then, the problem (1.2)-(H) supplemented with boundary 
conditions (1.5) and initial data satisfying (1.6) admits a weak solution {g,\i,ri\ on 
(0,00) X m the sense of Definition 2.2. In addition, 

i) the total fluid mass and particle mass given by 



Mg{t)^ I g{t,-) dx and M^{t) = / 7j{t,-)dx, 
Jn Jn 

respectively, are constants of motion. 
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ii) the density satisfies the higher integrability result 

g e L'<+^{{0,T) X n), for any T > 0, 
where 9 = niin{|7 — 1, |;}. 



We can now completely characterize the large time behavior of free-energy solu- 
tions to (1.2)-(1.7). 



Theorem 2.2 (Large-time Asymptotics) Let us assume that (fi, $) satisfy the 
confinement hypotheses (HC) . Then, for any free-energy solution (g, u, rj) of the prob- 
lem (1.2)-(1.4), in the sense of Definition 2.2, there exist universal stationary states 
gs{x), r]s(x), such that 

git) — > gs strongly in L'^{fl), 



esssup / p(r)|u(r)|^ 0, 

rjit) rjs strongly in LP^(fl) forp2 > 1, 



as t — > oo, where (77s, f?s) fl?^e characterized as the unigue free-energy solution of the 
stationary state problem: 



given by the formulas 



gs dx ^ go dx, rjs dx ^ rjo dx. (2.7) 
n Jn Jn Jn 



7-1 
a7 



h/31> + C,)] + 



r]s = C,,exp(-$), 



where and Cg are uniguely given by the initial masses due to (2.7). 



3 Global-in-time existence 

This section is devoted to the proof of the existence result (Theorem 2.1) which will be 
achieved by patching local-in-time solutions with the aid of suitable global estimates. 

Lemma 3.1 The conclusion of Theorem 2.1 holds true on any time-space cylinder 
[0,T) X ri, where T > is any given finite time. 

Taking Lemma 3.1 for granted, a weak solution verifying Theorem 2.1 can be 
constructed as follows. 
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Proof of Theorem 2.1: Fix any To > 0. From Lemma 3.1, we have the existence 
of a weak solution (pi, Ui, 771) on [0, 2Tq). To proceed, we wih need the fmiction 



fl, te[0,To-e], 
CW= <^dist(t,To), (To-e, To), 
[0, otherwise. 

By setting (pQ'^ as a test-function in the continuity formulation (2.3), we obtain 



/ / B{gi)ipe ^ dxdt 
J To-e Jn 

I " In (^'^(^1)^*'^ + CB{gi)ui ■ V,(/? - C&(f?i)div.Ui^) dx dt 



'0 Jn 

b{qoMo,-)C{o) dx, y^eCT{[o,T]xTl). 

'n 

Sending e — > in this formulation yields 



J J (^B{gi)dtip + B{gi)ui ■ Vx^fi - b{gi)divxUi(p^ dx dt 

B{g{To,-)MTo,-) dxdt^ f B{goMOr) dx, ^g e C^{[0,T] xTl). 



Hence, (gi, Ui) is a renormalized solution of the continuity equation on the closed in- 
terval [0, To] with additional boundary terms at t = Tq. By applying similar arguments 
to the momentum equation (2.4) and particle density equation (2.5), we conclude that 
(gi,ui,77i) is a weak solution on [0,To] (with additional boundary terms at t = To). 

By the various uniform in time bounds (cf. (2.6)), {gi,Ui,r]i) at time Tq has 
sufficient intcgrability to serve as initial data for a new solution (g2, U2) ^2) defined on 
[To, 2To] X ri. The energy of the new solution is less than the initial energy E{go, uq, ?/o) 
and the triple {g, u, r/) given by 

{g,u,i-i){t) = (ft,u,,77,), t e ((i - l)To,iTo], 

is a weak solution on [0, 2To] x fl. A solution for all times is readily obtained by 
iterating this process. □ 

In the remaining parts of this section we prove Lemma 3.1. The overall strategy 
can be outlined as follows: First, we prove that there exists a sequence of approximate 
solutions {g5,h,'ns.h,'r]s.h}h>o (see the ensuing subsection for details). Then, by first 
sending ft, — > and subsequently 5 — >■ in this sequence, we prove that in the limit 
we obtain a weak solution to the system (1.2)-(1.4) in the sense of Definition 2.2. 

3.1 Approximation scheme 

The approximation scheme is realized in two layers. In the first layer we add a 
regularization term to the pressure. 

Definition 3.1 (Artificial pressure approximation) For a given 5 > 0, we say 
that the triple {gs,us,r]s} is a weak solution to the artificial pressure approximation 
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scheme provided that {gijU^jT;^} is a weak solution in the sense of Definition 2.2 hut 
with the pressure psiffs) given by 

psigs) =p{gs) + Sg% 

and initial data {gsfl,mQ^s,r]sfi) ~ {go,mo,riQ) -k ks, with ks denoting the standard 
smoothing kernel. 

The second layer of approximation is a discretization of the equations in time. To 
define this layer of approximation, we shall make use of the space 

w(n) = L\n)r\L^{n) x w^ '^in) x w^■^n)nL^{n)r\L^{n). 

Definition 3.2 (Time discretization scheme) Let S > be fixed. Given a time 
step h > 0, we discretize the time interval [0, T] in terms of the points t^ = kh, 
k = 0, . . . , M , where we assume that Mh = T . Now, we sequentially determine 
functions 

{e^;,„<„%'Je fc = i,...,M, 

such that: 

• The time discretized continuity equation, 

d,^[e,%] + div,(^.i,uL) = o, (3.1) 

holds in the sense of distributions on Q.. 

• The time discretized momentum equation with artificial pressure, 

d-t [gih-^ih] + d^^Agih-^ih ® Kh) - t^^K.h - Av.div,^;, 
+ V. (psigih) + vlh) = ~{Pgih + %\)v.$, 

holds in the sense of distributions on fi. 

• The time discretized particle density equation, 

d^[r,i^] + div, {rilM,h - V.$)) - A%',. - 0, (3.3) 

holds in the sense of distributions on Q,. 

In the above equations, [0'^] — denotes implicit time stepping. 

For each fixed ft- > 0, the time discretized solution {gi^',^, ^, is extended 

to the whole of (0,r) x by setting 

{gs,h,nsj„ri&.h){t)^{gih.Vilh,rilh). t e {t''-\t^], k = l,...,M. (3.4) 
In addition, we set the initial data 

igs,hiO)Ags,o'a.s,h)iO),Vsji{0)) = (£'a\o,toa\o,%o)- 

Next two subsections treat the existence of a well-defined approximation scheme 
in both relevant cases 51 bounded or not under the confinement conditions (HC). 
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3.2 Well-defined Approximations: Q bounded 

To prove existence of a solution to the time discretized approximate equations (3.1)— 
(3.3) we will utilize a fixed point argument. For this purpose we define the operator 

Th [■] : w{n) ^ wm, 

as the solution (g, u, rj) = Th [''/', z, C] to the system of equations: 

+div,.(gu) = 0, (3.5) 

gu — Q^^^VL^^^ 

h div3,(£iu(K) u) - /iAu - AVi^div^u 

= -V:,(w(f?) + C)-(/3f^ + C)V.<i>, 



and 



"L^tl + div, (,7(z - V,*)) - Ar; = 0, (3.7) 



in the sense of distributions on VI, where [g^^^ ,u^^^ ^if'^^) G W{^) is given. 

Observe that for a fixed /i > 0, (5 > 0, and fc = 1, . . . , M, any fixed point of the 
operator Th [■] will be a distributional solution to the time discretized approximate 
equations (3.1)-(3.3). i.e the fixed point 

{BXh,^ih. vs,h) = T h [gih.'^ih, vs,h] 

is precisely the desired solution at each time step. 
3.2.1 is well-defined 

The existence of functions {g,u,i]) € yV(r2) satisfying (p, u, 77) = Th [Vs^i, C] follows 
from Lemmas 3.2 and 3.4 below. 

Lemma 3.2 ([18], Theorem 6.1) Let h > be fixed arbitrary and assume that 
ip'^'-i € L''{n)nL^in), j^m''-^ e L^n)nLP{n) for somep> 1, C e W^'^iil), and 
Va;$ G LP{n), for all p. There exists a weak solution {g,u) € LT(17)nL'5(rj) x '^(rj) 
of (3.5) -(3.6) satisfying 

ps{g) € L^n), (3.8) 



V^curl^u, V:, l^div^u - j^psig)j , G i''(A'), (3.9) 

where 9 = jf > f ^'^'^ A' C is any compact subset. 

Notation.- Throughout the paper we use overbars to denote weak limits, just to 
illustrate this notational use, we will often use convexity arguments based on the 
following classical lemma [11]: 

Lemma 3.3 Let O be a bounded open subset 0/ with M > 1. Suppose 17: M — > 
(—00,00] is a lower semicontinuous convex function and {fn}„>i is a sequence of 
functions on O for which Vn v in L^{0), g{v„) € L^{0) for each n, g(vn) 
g(v) in L^{0). Then g(v) < g{v) a.e. on O, g(v) G L^(0), and Jpgiv) dy < 
liminf„_i.oo /q <7(fri) dy. If, in addition, g is strictly convex on an open interval 
(a, &) C M and g{v) — g{v) a.e. on O, then, passing to a subsequence if necessary, 
Vn{y) v{y) for a.e. y G {y & O \ v{y) € (a, b)}. 



11 



Lemma 3.4 Assume that fl is bounded. Let if'^ e i^(f2) n (^2) n {77''"^ > 0}, 
and z G Wj'^(r2) be given functions. Then, for each fixed h > 0, there exists a non- 
negative function 77 g W^'^^il) n L^{U) satisfying (3.7) in the sense of distributions 
on CI. Moreover, 

[ /7-i[77log 77 + 77$] + |2V,V^+V^V,$p dx<C, (3.10) 

where the constant C > depends on ||z|| wi-2(o), ||77'^^"'"||l3(j7), and ||<l'||n/i,°=(n). 

Proof: For each e > 0, we let z' = z * k^, where * is the convolution product and 
is the standard smoothing kernel. From [21, Proposition 4.29], wc can assert 
the existence of a unique weak solution rj'^ G W'^''^{n), rf > 0, to the linear elliptic 
equation 

^^—^ + div:,, (r;'(z" - \7^<P)) - A77' = in n, (3.11) 

satisfying the boundary condition 

77'^Va;$ • ly + "^xlf -v^Q. 

By integrating over Q, (using the boundary condition), we observe that 

T]' dx= I T^-^ dx < C, 
n 



and hence 77*^ g L^{V,) independently of e. 
Let the sequence {Bi}]^^ be given by 



Bi{y) 



logy, y>l^, 
logri, y<r\ 



Moreover, let = {x £ U : r]%x) > l'^} and nf ^ n \ O;. Multiply (3.11) with 
Bi{r]'^) and integrate by parts to obtain 

/ ft."^ [77^ log77'^ - 77''~Mog77'^] dx - /i"^ [77*^ log/ - Ty'^^MogZ] dx 

+ / V.^$ • V:r?f + — |V:r??T " z'V^^rf dx = 0. 

Applying the identity xrj'^]'^ = 4|Va;y^p and reordering terms in (3.12) 

/ h-^T]" \ogr]^ + 4\WxV¥\'^ dx 
Jni 

= [ h-^ [t]' log I- ri''-^ log I] dx+ [ h-^r^^-^logT^' dx 
Jni Jill 

- / 2V^V:,$ • V^0^ - 2V^z^ • VxV¥ dx 
Jn, 



< I h-^\0glr dx+ I h-^Tj'^-^T]' dx 

Jni Jn, 

<h-Hoglr\n^\+ f h-\f-'ij' dx + 2\\V^\\l,^,,^^ 
Jn, 
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where the last inequality is the Cauchy inequality. The last term is bounded since 
'q^,r]^~'^ e L'^{Vl) and $ G M^^'°°(ri), independently of e and I. Using the Sobolev 
embedding W^'^(O) C L^(ri) together with jy''"^ e i^(f^), we achieve 

I /i-i77Mog77^ + 2|V,V^|2 da; 

(3.13) 

<C+ / ry^/i^V-' + dx<C{l + \W\\^. ). 

Applying the Holder inequality, Sobolev embedding, and Young's inequality (with 
epsilon), respectively, we bound the last term as follows 

Inserting this expression in (3.13) and fixing /? small, we gather 

/ h-^Tj'\og7j' + \\/,^f dx<C, (3.14) 

where the constant C is independent of both e and I. Since is bounded, it follows 
that 

independently of e. 

Since W^'^{^1) is compactly embedded in ^^(ri) for all p < 3, we can conclude 
that, passing to a subsequence if necessary, 

ri'^^T] in W^'^n), ri'^^T] in LP(17), p < 3, 

as e — )- 0. Consequently, we can take the limit e ^ in (3.11) to conclude that 77 
satisfies (3.7) in the sense of distributions on J7. 

Next, since 77^ 1] in LP{n), p < 3, we can also conclude that 

Vv^^y/7j in L^iil), q<6, \7a:VTf -^"^xVv in L'^{^), 
as e — > 0. Using this, we can send e — > in (3.14) to obtain 

/ h^^f] log T] + \Wxy/v? dx < C. 
Jn 

By weak lower semi-continuity || Va;y^||^2(Q) > !I^sv^IIl2(j2)- Thus, 

Jn 

Then, using the identity 

we are led to the conclusion 

/ /i-i [?7log?7 + ?7$] + |2V^^Af + V7V:,$p dx 
Jn 



<C+ [ /^-^/$ + 4V,V^- (y^V,$) + ?7|V,$|2 
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which is (3.10) and the proof is complete. 



□ 



Remark 3.1 Since V xs/rj € 



(f2), the set of particle vacuum regions have measure 



zero; \{x G il; = 0}| = 0. 
3.2.2 Th admits a fixed point 

Wc now prove that Th [■] admits a fixed point and consequently that the time dis- 
cretization scheme in Definition 3.2 is well defined. The key observation made is that 
the bound on the pressure enables us to obtain an energy equality. This equality 
in turn yields compactness of the operator Th [•]• 

Lemma 3.5 Assume the case of hounded domain Vl. Let {g'"'^^, jti'"'^^, rf^^} G i^(il) x 
Z/5 (n) n L^ip.) X L"'{^),p > 1, be given functions. Then, for each fixed h > there 
exists a fixed point {g, u, ri} G W{fl) for the operator Th [•]; i-e 



As a consequence, the time discretization scheme given by Definition 3.2 is well- 
defined for bounded domains Q. 

Proof: We will prove the existence of a fixed point by verifying the postulates of 
the Schauder corollary to the Schaefer fixed point theorem [9]; // the operator Th [■] 
is continuous, compact, and the set {s G [0, l],a; G yV(fi) : x = Th [sx]} is uniformly 
bounded, then the operator Th [■] admits a fixed point. 

First, we observe that the operator Th [•] is clearly bounded and continuous. 

Next, we prove that the operator Th [•] is compact. For this purpose, let {0, z„, 77„}^ 
be a sequence such that {zn,r]n) Gt M^o'^(r2) x W^'^{Q) for all n = l,...,oo, and 
construct a sequence {gn, u„, 77„}^j^ by setting 



Then, from the previous lemmas we know that (gn, u„, ?7„) G yV{^l) independently of 
n. Hence, we have the existence of functions {g,u,ri) G W{fl) such that, by passing 
along a subsequence if necessary. 



Moreover, by compact Sobolev embedding, we clearly have the existence of (0, z, 77) G 
yV(r2) such that — > C 9-6 in ^ and z„ — > z a.e in fJ. 

Now, we claim that in fact (p„,u„,77„) — > {g,u,i]) in W{^l), where (g, u, 77) = 
Th [0,z, C], and consequently that the operator Th [■] is compact. 

To prove this claim, we first note that compactness of 77„ in W^^' 2 (fi) is immediate 
from the linearity of (3.3). That is, since r]n G W^'^ {fl) we have by Sobolev embedding 
that rjn ^ rj a.e in Vl and thus by setting log 77— log rjn (cf. Remark 3.1) as test-function 
in (3.7) one discovers 



{g, u, r/) = Th [(?, u, 77] . 



(Pn, u„, ??„) = Th [0, z„, Cn]i n = 1, ... ,00. 



(e„,u„,77„) (e, u,77), in yv(ri). 





(n). Moreover, in the limit we have 



(3.15) 
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in the sense of distributions on O. 

We continue proving compactness of the operator Th [•] by first proving strong 
convergence of the density Qn ^ Q a-c in f7. Compactness of the velocity u„ in 
VFq '^(£1) will then follow from an energy equality. In order to prove strong convergence 
of the density we will need weak sequential continuity of the effective viscous flux. 
That is, 

lim / [(A + ^)div^u„ - ps^Qn)] i^Qn dx 

^ (3.16) 

(A + ^t)divj,u-p5(£() i^gdx, yip € Cf [VL). 

Here, (3.16) is an immediate consequence of (3.9) and compact embedding of Sobolev 
spaces. 

Before we proceed to prove strong convergence of g„, we first note that there is 
no problem with taking the limit in (3.1) to obtain in the limit 

^Q + AWx{\lQ)^^Qk-l, (3-17) 

in the sense of distributions on il. Hence, since in particular g E L^{il) and u e 
Wq''^{Q) we can conclude that, for any B G C[0, oo) n C^(0, oo), 

^gB'ig) + div,(B(g)u) + {(gB'ig) - B(f?))div,u) = ^gk-iB'{g), (3.18) 

in the sense of distributions on fl. Then, by setting B{z) = zlogz we obtain the 
identity 



gn log gn ~ g\ogg+ gn - g dx 

= / £»div^,u - £)„div^u„ + gk-i{log g^ - log g) dx. 
Jn 

Thus, by taking the limit n — > c», we see that 



g log g — g log g dx 

. . (3.19) 

= lim lim / i/jig div^u - ipigndiv dx + / gk-i{\og g - ^og g) dx, 

where i>i G C^{Q)n{ipi > 0} is such that ■0/ = 1 on the set {x G fi; dist(a;, 9ri) > j}- 
By an application of (3.16) we see that for each 1^1... oo. 



lim / tpigdiv^u- ipigndiv^Un dx = lim / 'il)i{ps{gn)Q ~ Ps{Qn)Qn) dx<0, 

where the last inequality follows from the convexity oips{gn)- Similarly, the concavity 
of log gn yields 

/ gk-i{logg -\ogg) dx <Q. 
Jn 

However, then (3.19) together with the convexity of zlogz gives 



0< / g\ogg — gloggdx < 0. 
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which immediately yields p„ — > g a.e in f2 due to Lemma 3.3. 

We are now in a position to prove that u„ — >■ u in Wq''^{Vl). For this purpose, we 
first set u„ as test-function in (3.2) to obtain the identity 

hm / — 'r li\V^\in\ +\\diY^Vin\ dx 

n-s-oo Itl li 

= lim / (psign) + Cn) diVx^n ~ (Cn + /3e'n)Va;^> • u„ dx (3.20) 

= / (p5(^?) + C)div^u-(C + ,3^?)V^$•uc^a;. 

Note that in addition to the strong convergence of ?7„ and g„ we really need (3.8) to 
conclude this. 

Now, since gi„ Q, Cn ~^ Ci ^-nd u„ — > u a.e in there is no problem with passing 
to the limit in (3.2) to obtain 

^ ^ diva;(£iu u) — /lAu — AVjjdivjjU 

= -v,(w(^^) + C)-(/5^' + C)v.$, 

in the sense of distributions on Vl. Hence, in view of (3.15) and (3.17), we can conclude 
that (g, u, 77) = Th [0, z, C]. Moreover, by using u as a test-function for this equation 
we obtain the identity 



't ^Ji\V ^-mY + A|div^u|- 



dx 



= / (p5(^?) + C)div^u-(C + ^g)V,.$•uc^a; 

= / ;r7 1 h /i VxUp! + A div^u^dcc, 

Jn 2/i h 

where the last equality is (3.20). This can only happen if Va;U„ — > WxU. Thus, we 
can conclude that T/i [■] is a compact operator. 

Finally, let s G [0, 1] be arbitrary and assume that there exists a triple (p, u, 77) e 
yV(ri) such that (g, u, rj) — Th [sg, su, s-q]. Then, by setting u as test-function in (3.2) 
we get the identity 



f gu- + gk-iu'' unik-i , ,2 , ,2, 

/ TTT 7 h M V^up + A diva^upda; 

Jn 2/i h 

= / (Psig) + S77)div^u - {sri + /3e)V^$ • u dx. 
Jn 



(3.21) 



Using B{z) = 115(2) := Tpjp{g) + f as the rcnormalization function in (3.18) we 
also have the identity 

- [ psig)divxu dx ^ [ Wgig)ig - gk-i) dx. (3.22) 
Jn " Jn 

Similarly, using /3$ as test-function in (3.5) gives 

^ g- gfc-i ^ - /^g u • V^<^>dx ^ 0. 
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(3.23) 



Using $ + log 77 as a test-function in (3.7) (cf. Lemma 3.4) gives 

= / - — ($ + log ?/) - sr/u • V^:* + S77div^u 
Jn h 

+ |2V:,V^+0yV^$p dx. 
By combining (3.22)-(3.23), wc obtain the identity 

/ Ps{(})divxU + sTjdivxU — {st] + (3g)S/x^ ■ u dx 
Jn 

f / (g- gfc-i) 1] - rjk-i 
^Jn '^'^ h ^($ + log^) 

Then, inserting this into (3.21) and reordering terms gives 

gu^ + Qk-iu^ U77ifc_i g-Qk-i n^. , V-Vk-i,^ , , ^ , 
2^^ — + ("^-^^^ + + + ^^S'^)'^^ 



+ / ^i|Vj;Up + A|div^up da; 
Jo 

/ \2\/,^+^V,^\^dx<0. 
Jn 



(3.24) 

Consequently, we can conclude that 

where the constant C depends on the data tjk-i, Qk-i, ^k-i, and $, together with 
h. However, C does not depend on the parameter s. 

We can now conclude the proof since we have proved that the operator Th [•] is 
bounded, continuous, compact, and that the set {s g [0, l],x S W{il) : x ~Th [sx]} 
is uniformly bounded. □ 



3.3 Well-defined Approximations: ^2 unbounded 

At this point we have proved that the approximation scheme is well-defined on 
bounded domains. Given an unbounded domain and an external potential <i> satis- 
fying the assumptions (HC), we can always find an increasing sequence of domains 
Qr, with r > such that fir are bounded and (flr,^) satisfies (HC) approximating 
fl in the sense Ur>Q^r = Using the previous subsection, for any r > 0, there is 
a solution on fir- In this subsection, we prove that we can send r — >■ 00 to obtain a 
solution in fl. The following lemmas will be of use in the sequel. 

3.3.1 Consequences of confinement: unbounded 

We show in this part how to control the negative contribution of the physical entropy 
?7log7y in the free-energy bounds for unbounded domains 51. Here, the confinement 
conditions (HC) on {fl, $) are crucial. Most of these lemmas can be seen in [10] but 
we include them here for the sake of completeness. We first start by a classical lemma 
in kinetic theory. 
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Lemma 3.6 Assume that (rj,$) satisfy the hypotheses (HC). For any density rj G 

r]{x) log_77(x) dx<- [ <^{x)t]{x) dx + - [ e-*^^)/^ da; . 
Proof: Let fj := ?/X{,,<i} and M = fj{x) dx < 7]{x) dx — M. Then 
fj{x) (\ogfj{x) + l-^(x)] dx::^ I [U {x) log U (x)] fi dx - M log Z 



where U :— fj,{x) = e with Z ^ J^e dx. The Jensen inequahty 

yields 



[U{x) logU{x)]ndx > / U{x)ndx\ log / U{x)ndx] =Af logM 



and 



ri{x) log_ ri{x) dx ~ j f]{x) logfj{x) dx > MlogM — MlogZ — [ ^{x) fj{x) dx 

Z 1 /" 
> / r;(a;) , 



n 



e 2 

from which the desired claim follows. □ 



We can immediately use this previous lemma to conclude the following conse- 
quence. 

Corollary 3.1 Assume that (f^, $) satisfy the hypotheses (HC). For any density 

V e LUn), if 



T]{x) logri{x) dx + / ^{x)r]{x) dx < C , 
! Jn 

then 77 log 77 G L^{fl) and there exists D > depending on C and <i> such that 

r]{x) log_|_ ij{x) dx < D and / ^{x)r]{x) dx < D . 
n Jn 

Finally, the above estimates can be used to control the mass of the densities rj 
outside a large ball to avoid loss of mass at infinity. 

Lemma 3.7 Given any domain fl such that e^* £ L^{il) and any density rj G 
L\{n), then 



ri{x) log ri{x) dx + / ^{x)ri{x)dx> / 77(0;) da; log 



In '?(^) 

fn Jn Jn v/^e-^^da; 

As a consequence, if e^^ E L^{^1) and 

r]{x) logri{x) dx + ^(x)ri{x) dx < C , 
n Jn 

then, for any e > there exists R > depending on C and $ only such that 

r]{x) dx < e . 

nn(R3/_B(o,i?)) 
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Proof: A direct use of Jensen's inequality shows the first inequality by using the 
convexity of x i— > a; log a; . To show the second claim we start by applying the first 
inequality to the domain fi^ := n (R^ /B{0, R)) from which we obtain 

rjix) dx log < D (3.25) 



for some D > 0, where Lemma 3.6 and Corollary 3.1 were used. Now, we argue by 
contradiction, if the second claim were not true, we would have 

3e > Vi?o >03R> Ro such that / r){x) dx > eq. 
Since e~* £ L^{n), we can always assume that Rq is large such that 
e-*(^)da;< / er'^^'^Ux < eo < f ri{x) dx 



and thus due to (3.25), 

r,{x)dx< I e-*(^) da; e-°/^° < / e-'^'^^'Ux e^/'° . 



This leads to a contradiction since the right-hand side can be made arbitrarily small 
by taking Rq large enough. □ 



3.3.2 The approximation scheme is well-defined 



Notation. - In what follows, we will often obtain a priori estimates for a sequence 
{''^n}n>i that we write as "w„ X" for some functional space X. What this really 
means is that we have a bound on that is independent of n. 

Lemma 3.8 Set(fl,^) satisfying (HC) withQ unbounded and let {g'^~^ ,u'^~^ ,r]^~^) € 
W(ri) be given data. The time discretization scheme (3.1)-(3.3) admits a distribu- 
tional solution in the sense of Definition 3.2. 

Proof: Let {i7r}r>o be an increasing sequence of domains such that Ur>o^^r = ^ 
and such that, for each fixed r, Vlr is bounded and ($7^,$) satisfies (HC). From 
the results of the previous subsection, we have the existence of a triplet (prjUn'yr) 
satisfying the time discretized equations (3.1)-(3.3) in the sense of distributions on 
r^r. Consequently, we can define a family of such solutions: 

{£lr,Ur,77r}_R<r<oo, (gr , Ur , 77r ) € yV(rir ) , for Cach fixed T S (i?, Oo) , 

where R is fixed according to the requirements on the potential (see (2.2)). For this 
construction, (3.24) yields 



QrK 



2h 



f (n',(^,,) + /?$) 



julVajUrp + A|diva;Ur 



h 



log?7,.)da; 



< 



< eh' 



(3.26) 



h ^rjk-i logry^ dx 



1 



\m-i\\mn) + ll^fc-il 
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In addition, integrating (3.3) over fir, 

with constant C independent of r. An interpolation inequality and the Sobolev- 
Poincare inequality allow us to conclude 

13 3 

llVrh^n^) < \h-\\li^nr)\Ml^n^-, < C\\^xVv^\\l2^n^y (3-27) 

where the constant C is independent of r. By applying (3.27) and the Sobolev- 
Poincare inequality to (3.26), we conclude 

i '-^=^^^i^+'ii^'s{9r)+m + ^-^{'^ + \ogr^r)dx (3.28) 

+ / (//-e)|V^u^p + A|div^u^|2 + |2V^7^+V^V:,$p-e|V^V^pda;<C, 

where C is independent of r. 

Next, we observe that (3.28) yields 

/ 4|V,.0^p + 2V^?7,.V^$ + 77r|V:,4>p = / |2V^V^+ V^V^$p d.T < C. 
Reordering terms and integrating by parts 

/ 4|V:r7r^p +77r|V:r$p dxdi < C+ / ?7r|A$| da;. 
Then, (2.2) and (3.28) gives 

< ||A$||ioo(n«) / r]rdx + C 77^* dx < C, 



with C independent of r. 

Setting (3.29) into (3.28), fixing e smaU, and applying Corollary 3.1 gives 

i + g'"'"' + 1 (n- (g.) + /3ci>) + ^(ci> + log^ ^r)dx 

+ / ^|Va,Ur|^ + A|div3;Ur|^ + iVa^V^I^ + 'yrlVjj^l^dcc < C. 
Jn,. 

Since Ilv^lli2(a ) ~ In Vr dx < C, the previous estimate allow us to conclude 

r]r eb L^i^r)nW^'i{nr), (3.30) 

Next, let (0, V, ?/;) € [Cjj"(51)]^ be arbitrary and fix a number r G (i?,oo) large 
such that 

supp((^, V, ■(/)) C fir- 
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Then, from the previous bounds, we have the existence of functions {g, u, i]) G yV(ri) 
such that, along a subsequence. 



(3.31) 



gr g, in L^(fip), ^ u, in W^''^{rtf), 

Moreover, by compact Sobolev embedding 

rjr'^'n, in Lfoc(f^f), ^ u, in ^[^^(f^f), p < 6. (3.32) 

Using (3.31), (3.32), there is no problem with sending r ^ in 

dt[gr]'t> ^ grUr'^x4> dx — 0, 
/ dflrir]^^ + {Ur - Va;$)Vj;?/' - Va;7?rV:rV dxdt = 0, 

to conclude that (gi, u,7y) solves (3.1)-(3.3) in the sense of distributions on 57 (Recall 
that (0, ip) was chosen arbitrary). Similarly, we can pass to the limit in the momentum 
equation (3.2) to obtain 

/ dj'[giu]v — pu (g) u : V^jV + /iVj^uVj^v + Adiv^jUdiv^-v dx 

+ / {[3 g + x^'v dx ^ lim / ps{gr)divx\' dx. 

Hence, it remains to prove that ps{gr) ^ Ps{q) as r — > oo. This problem is also the 
core problem encountered in the existence analysis of Lions [18]. Due to the regu- 
larity properties of rjr, the presence of the rj variable does not impose any additional 
difficulties as compared with [18]. However, some minor modifications are needed to 
incorporate the unbounded potential The needed modifications for the present 
stationary case are almost identical to those of the non-stationary case. Thus, we do 
not give the arguments here but refer the reader to the more general arguments given 
in Section 4.5.2. □ 



3.4 There exists an artificial pressure solution {h — )• 0) 

In the previous subsection we proved that, for each h > Q and 5 > Q, we can con- 
struct functions {gs.h, ^s.h^ Vs.h) according to Definition 3.2 and (3.4). In this section, 
we prove that the corresponding sequence {{gs,h,Vis^h,Vs,h)}h>o, with S > fixed, 
converges as the time step /i to an artificial pressure solution in the sense of 
Definition 3.1. 

Lemma 3.9 [2!3, Corollary 4, p. 85]. Let X C B C Y, be Banach spaces with X C B 
compactly. Then, for 1 < p < oo, {w : u G LP{0, T; X), Vt G L^{0, T; Y)} is compactly 
embedded in LP(0, T; B). 

The following lemma is a variation of a result due to P. L. Lions (see [17]). 



21 



Lemma 3.10 For a given T > 0, divide the time interval (0,T) into M points such 
that (0,T) = U^£]^ tfe], where tk = hk and we assume that Mh = T. Let 
{Ai/oO' {9h}'h>o ^^"^ sequences such that: 

• {fh}h>o,{9h}h>a converge weakly to f,g respectively in LP^{0,T]L'^^{fl)), 
iP^(0,T;L92(rj)) where 1 <pi,qi < oo, 

1 11 1 _^ 
Pi P2 qi 92 

• the mapping t — >■ gh{t,x) is constant on each interval {tk-i,tk], k — 1,. . . ,M. 

• the discrete time derivative satisfies, 

9^^^^-)-9H{t-h,x) g^^i(o^r.^-M(^)). 

{fh]h>o satisfies: ||/« - /„(■ +C, OIIlp2(o,T;L<!2 (o)) 0, as \£,\ -> O,uniformly in 



n 



Then, gtfh ^ gf 'in the sense of distributions on (0,T) x fl. 
3.4.1 Energy estimates 

Let (5 > be fixed and let {gs.h,'*i5,h,ils.h}h>o be a sequence of time discrctized 
solutions constructed according to Definition 3.2 and (3.4). Since (3.24) holds for 
every k, we can sum this equality over all A: = 1, . . . , m, for any m € [1, M], to obtain 
the energy equality 

+ / nl^xUsji]'^ + \\div^us,h\'^ dxdt 

Jo Jn 

+ f f |2V^7^+7??ZIV^$p dxdt 

Jq Jn 

rn 

+ V / %-ilog 



dx 



fe=i • 

where t"^ = mh G (0,T) and the energy E{-,-,-) is given by (2.6). By convexity 
(7 — l)a'>' + f — 'fa'^~^b > for all a, 6 > 0, 7 > 1. By concavity of z 1— )■ log z, 

ijk-i log f ^) > / 77'"' -v" dx^ 0. 
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Due to the confinement conditions, this, Corohary 3.1, and (3.33), allow us to conclude 
the bounds: 



gs,h ef,£°°(0,r;iT(i])nL6(r!)), 
m,hiogT^s,hebL^io,T;L\n)). 



(3.34) 



By similar arguments as those leading to (3.30), we conclude 

Vs.h €6 L\0, T- L\n)) n L\Q, T; W^'^n)), (3.35) 

independent of both h and 5. 

Utilizing the above h independent bounds, it is a simple exercise to obtain from 
the system (3.1)-(3.3) the weak time difference bounds 

d't'[gs,h] ebi°°(0,r;M^-i'i(f])), 

dlim.h] L\0,T;W~'^Hn)), (3.36) 



3.4.2 Convergence 

In view of (3.34) and (3.35), we have the existence of functions {gs, u^, rjs) such that, 
along a subsequence as /i — > 0, 

g5,h^g5 mL^{0,T;L"'{n)), 

m,h^m mL^{0,T;L\n)). 
By virtue of (3.36), Lemma 3.10 can be applied to yield 

gs,h^sji gsus, g5MUs,h ® u^^/i ^ gsus ® u^, (3.37) 

in the sense of distributions on (0, T) x as — > 0. Here, gs,hUs,h ® us^h gs^s ® us 
follows from setting gh = gs,hU.5^h and fh = Uh^s in Lemma 3.10, where gh ^ 9 = gs^s 
from the second last result in (3.37). 

Next, since 775,,, ^b L\0,T; L\n))nL^{0,T;WHn)), dtVSM €b L\0,T;W-^^\n)) 
and is compactly embedded in L^, we can apply Lemma 3.9 to conclude that 

ris.h^r^s inL2-^(0,r;iL(^^)), V^I ^ Vv^ in LH0,T; Ll^{n)), (3.38) 

Now, by a straight forward application of the Holder inequality, we deduce \'x'r]s,h Gfo 
L^{0,T; L^{U)) n L^(0, T; (ri)). The standard interpolation inequality provides an 
estimate of the form WxVS,}! &b L°'^{0,T;L°'^{fl)). Thus, in particular 

'^xVs.h ^ '^x'nsi (3.39) 
in the sense of distributions on (0,r) x fl. 
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Now, equipped with (3.37), (3.39), and the bounds (3.34), there is no problem 
with taking the Hmit ft, in (3.1) and (3.3) to discover that 

- / / gsiipt ~ usyx-ip) dxdt ^ / Qs,oijiO,x) dx, (3.40) 
Jo Jn Jn 

for all V e C^i[0, T) x fl)) and 

775 (V.'t + [us - V^<^)V:,ip) - Vxm'^xip dxdt ^ [ ?75,o'0(O, a;) dx (3.41) 



Jn 



for all e C^{[0,T) x Q)). 

Hence, the limiting functions satisfy both the continuity equation and particle 
equation in the sense of Definition 3.1. 

Similarly, we can go to the limit /i ^ in (3.2) to discover in the limit 

/ gsugvt + gsus ® : V^^v^ + r/^div^^v dxdt 
Jn 

= / / fi'VxUs'VxV + XdivxUsdivxV + {Pgs + ris)^^^^ dxdt 
Jo Jn 

— lim / / ps{gs h)^^^^ dxdt — / mov(0, x) c?x, 
'»^oJo Jn ' Jn 

for all V S C^([0,T) x i7). Thus, in order to conclude existence of an artificial 
pressure solution in the sense of Definition 3.1, we must prove that in fact 

lim / / ps{g5ji)divxV dxdt = / / ps{gs)divxV dxdt. (3.42) 
^^ojo Jn Jo Jn 

Lemma 3.11 Fix any 6 > and let {gs,hT'^s.hi'>]s,h}h>o be a sequence of time dis- 
cretized solutions constructed according to Definition 3.2 and (3.4). Then, there ex- 
ists a triple {gs,us,'t]s) such that as /i — > 0, gs,h gs in L°°(Q,T;L^{^1) n L^{^)), 
gs,h gs a.e in (0,T) x D,, us,h in L'^{0,T;Wq' (fl)), gsjiVLs.h gs^s in the 

sense of distributions on (0,T) x Vl, rjg^h Vs in and rjs,h m 

a.e in {0,T) x fl, where {gs,Us,ris) is a weak solution to the artificial pressure ap- 
proximation scheme in the sense of Definition 3.1. 

Proof: In view of the high integrability, and strong convergence properties, of 
?7(5,/t, (3.42) can be proved by the same arguments as those leading to Theorem 7.2 
in [18]. Some minor modifications are needed to treat the unbounded potential <&. 
The arguments needed are identical to those of Section 4.5.2 and will not be given 
here. From this and the previous results of this section we can conclude existence of 
an artificial pressure solution. It remains to prove the energy inequality (1.8). 

We start with the following calculation 

lim / / 4|V^7^|2 + 4V^775.,,V^$ + ?75,h|V^$|2 dzdt 
h^oJo Jn 

= f f rM|V,$|2 +4V,7mV,$ + 4|V,VwP dxdt (3.43) 
Jo Jn 

> I I |2V,V^+ dxdt, 
Jo Jn 
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where we have used Lemma 3.3 and (3.39). By taking the Umit ft, ^ in (3.33) (using 
convexity oi z i-^ z'^ and z i— >■ zlogz), we obtain 



E{Qs,us,ris){t) 

+ n\VxUsf + X\divooUsf dxdt 

Jo Jn 

+ 11 |2V^v^+ V^V:,$|^ dxdt 
Jo Jn 



(3.44) 



/o Jn 
< E{gsfi,usfi,r]s,o), 

for any t e (0,T). □ 

3.5 Vanishing artificial pressure limit (6 — )■ 0) 

In the previous subsection we proved that, for each fixed (5 > 0, there exists an 
artificial pressure solution in the sense of Definition 3.1. Throughout this section, we 
let {gs,us,T]s}gyQ be a sequence of such solutions. The aim is now to prove that this 
sequence converges as (S — > to a weak solution of the fluid-particle interaction model 
(1.2)-(1.4) in the sense of Definition 2.2. This will then conclude the proof of Lemma 
3.1. 

3.5.1 Energy bounds 

The energy inequality (1.8), together with Sobolev embedding and Corollary 3.1, 
allow us to conclude the following (5-independent bounds 

gslusl"" Eb L°°(0,r;ii(r!))nL2(o,r;L"i(f^)),mi > 1, 

us Gfc L2(0,T;<^2(r!)), 

ei-efcL°°(o,r;i^(r!)), (3.45) 

gsusebL^{0,T;L^{n)), 

miogm ebL^{o,T;L'{n)). 

By straight forward applications of the Holder inequality, using (3.45) and Sobolev 
embedding, we deduce the bounds 

gsus Gb L2(0,r;i™^(l])), ggug ug €b L^{Q,T-L'-{n)), 

where 

67 37 

"^2 = — , C2 



6+7 3+7 
From (3.35), we also have that 

i^s &b L^{0,T;L^in)) n L\0,T;WHn)) 

Using the above (5— independent estimates, we easily deduce the weak time control 
bounds 

dtVsebL'{0,T;W-'^Hn)), 
dtgs ebL^{0,T;W-''^{n)), 
dtigsus) eb L\0,T;W-'-\n)). 



25 



3.5.2 Convergence 

The bounds in the previous subsection asserts the existence of functions {g, u, rj) such 
that, passing to a subsequence if necessary, 

gs^gs mL^{0,T;L^{n)), 

r,s^m inL^O,T;L^n)). 

As in the previous subsection. Lemma 3.10 can be apphcd to conclude 

gsus gu, gsus u^- gu<Siu, (3.46) 

in the sense of distributions on (0,r) x fl. 

By similar arguments as those leading to (3.38), we deduce 

and 

^xVs ^ V,77, (3.47) 

in the sense of distributions on (0,T) x fl. 

Using (3.46), the bounds (3.45), (3.47), and strong convergence of the initial con- 
ditions, we can take the limit 5 — )■ in (3.40) and (3.41) to discover that 

- / / Qii't - WVxip) dxdt = / eoV'(0, x) dx, 
Jo Jn Jn 

for all ^ e C^{[0,T) x H)) and 

- / [ V{ipt + {u- V2;$)V:j-0) - V^r/Va;^' dxdt = f miPiO, dx 
Jo Jn Jn 

for all ■0 € C^([0,r) X fi)). Hence, the limiting functions satisfy both the continuity 
equation and particle equation in the sense of Definition 2.2. 

Similarly, wc can go to the limit i5 — >■ in (3.2) to discover in the limit 

i-T I- 

guvt + gu(E)u : Vx'^s + f?diva;V dxdt 



Jn 



= / /iVa;uVa;V + Adiv:i,udiv2:V + {f3 g + 77)V2,<i>v dxdt 
Jo Jn 

— lim / / p{gs)divx'v dxdt — / m.ov{0,x) dx, 
■5^0 Jo Jn Jn 

for all V e C^{[0,T) x il). Thus, in order to conclude existence of a weak solution 
in the sense of Definition 2.2, it remains to prove that 

lim / / p{gs)divx'v dxdt = / / p{g)diVxV dxdt, 
"5^0 Jo Jn Jo Jn 

together with the energy inequality (1.8). Consequently, we are faced with a similar 
situation as in the previous subsection. The main difference is that we now only 
have 7 > Since rjs enjoys both high integrability and compactness, the proof 
follows by a small extension of Feireisl's arguments in [12]. First, wc establish higher 
integrability of the density on the entire domain fi. 
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Lemma 3.12 Let {gi^, u^, ry^j^-^Q he a sequence of artificial pressure solutions in the 
sense of Definition 3.1. Then, there exists a constant c(T), independent of S, such 
that 



I I qT^ dxdt < c{T), 



where 8 = niin{|7 — 1, j}. 

Proof: Since 775 Gb L^{0,T; L^{n)) (1 L^O,T■,W^'^n)), the addition of r;^ in the 
equations does not impose any potential problems. Consequently, if is bounded, 
the proof follows by the classical arguments ([18, 12]). 

If il is unbounded, then Vx$ is no longer integrable and we cannot simply apply 
existing results. To prove the bound in this case, let be the inverse Laplacian 
realized using Fourier multipliers (see [18, 11] for details). For each fixed S > 0, let 
the test-function ^fs be given as 

By the requirements on 9, we have in particular £b L°°{0, T; L''{Q)), 

Thus, V5 Gh L°°(0, T; W^-\^)) n L°°(0, T; L°°(fJ)). 

Next, since (p5,U5) is a renormalized solution to the continuity equations, (2.3) 
with B{qs) = states 

dtgt = -div,(^?^u) -{9- l)gU^v^u, 

in the sense of distributions on (0, T) x 57. For notational convenience, we observe 
that 

\\dt^s\\Lp(,Q,T,Li{n)) = \\^xA^^dtgs\\Lp{a,T;Li{n)) 

< \\9Ws\\LP{0,T;L<i{n)) + \\Q^<i^^xU\\LP{0,T;L^{n)), 

for appropriate 1 < p, 9 < 00 and r* ~ q. 

Next, we apply as test function for the momentum equation to obtain 



Jn 



agj^^ dxdt = — 



/ / (£'5U5)9tV5 + gsus ® Us : S^x^s dxdt 
Jo Jn 

T I- 

/iVa^u^VjjVa + Adiva;U5div2;V5 dxdt 



Jn 

T 



/ / 'nsQs - {qsP + 'ns)^x'^s^s dxdt - I moVi(0, •) dx 
jQ Jn Jn 



:= h+h + h. 

To conclude it remains to boimd /i, I2, and I^, independently of 5. We start with 
the Ii term: 

\h \ 0sdtv + Qs^s «> U5 : V^,v^ dxdt 

jQ Jn 

< \\gsus\\ _2j_ \\gs^s\\ ^ 

L°=(0,T:L7 + i (O)) L1(0,T;Lt-i (n)) 

+ ll£'5U5||i2(o,T:L™2(n))C(r)||/div^u||i2(o,r;L'-(n)) 

+ WQS^S <X) Us\\L2i^o^T;L'=2{n))\\'^x'Vs\\^2(^o^T:L^'2{n})' 
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where 

67 * / V I ^ 
r = , r = mo , Cn = < s. 

77-6' ^ ' ' ^ 27-3 - 

Now, we estimate 

||£)^div^U5||i2(o,T;L'-(n)) < C'l|div^Ui-||i2(o,T;L2(a))||£'5llL~(o,T;L=(0)), 

and hence conclude that < C{T). Next, we easily deduce the bound I/2I < C{T), 
and it only remains to bound I3. 

\h\ < \\vs\\L^o.T;L^n))\\g''\\ ^^^^ 

+ + ??5)Vx$|Ili(o,t;L1(o)) I!v5||l~(o,T;L°°(o)) + l|mo||Li(o) l|va||L=o(n) 
< CiT) (1 + Wi^gs + Vs)'^ML^io,T:LHn))) ■ 

Using the energy estimate (3.44), Corollary 3.1, and the requirements (2.2) on the 
potential, we readily deduce 



sup / \(3gs+T]s\\y.r^\ dx 

E(0,T) Jo 

< ||V:i,$||L«^(B(0,fl)) sup / 

te(o,T) Je 



figs + r]s dx 

e(o,T) Jb{o,r) 

+ C sup / iPgs + ?/5)$ dxdt < C{T), 

te{o,T) Jn\B(o..R) 

which brings the proof to an end. 

□ 

The following lemma concludes the proof of Lemma 3.1. 

Lemma 3.13 Let u^-, r^^j^^Q be a sequence of artificial pressure solutions in the 
sense of Definition 3.1. Then, there exists a triple {g,u,r]) such that as S ^ 0, gs ^ g 
inL°°{0,T;L'-'{n)), gs g a.e m (0,r)xf], ug ^ u in L^{0,T;WQ^{n)), gsus gu 
in the sense of distributions on {0,T) x Q, r]s ^ rj in and rjs ^ rj 

a.e in (0,r) x f2, where {g,\i,vi) is a weak solution to the particle interaction model 
in the sense of Definition 2.2. 

Moreover, the total fluid mass and particle mass given by 



Mg{t) ^ / git,-) dx and M^(t) ^ / ri(t,-)dx, 
Jn Jn 

respectively, are constants of motion. 

Proof: 1. Since rjs — > ?? a.e in (0,T) x fi, the presence of the rjs variable does not 
impose any extra difficulties as compared with the corresponding situation for the 
barotropic compressible Navier-Stokcs equations. Consequently, strong convergence 
of the density and existence of a solution follows as in [12, Section 4.2]. In the case 
Q bounded, this concludes the proof. For an unbounded domain, the potential is 
non-intcgrable and it remains to prove the energy estimate (1.8). 
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Since the mappings z ^-^ and z zlogz are convex, we can apply Lemma 3.3 
to conclude that, for a.e t £ (0,T), 



E{g,u,r]) < lim E{gs, Us, rjs). 

5— >0 



(3.48) 



Next, similar to (3.43), we deduce. 



/ / |2V,V^+ V^V,$|' dxdt 

JQ J K 

I I xy/m? + xvs^ + dxdt 

Jo JK 

I [ \2\/xy^+ ^/vs'^x'i'l^ dxdt 
Jo Jn 



< lim 

(5-s-O 



(3.49) 



= lim 

■5^0 10 Jn 

< E{go,Uo,r]o). 

for any compact K. 

Finally, using (3.48) and (3.49), we can send (5 in (3.44) to obtain 

E{Q,U,7]){t) 



/ / Ai|Vj;Up + A|div:rUp dxdt 
Jo Jn 

f f ry|V,.$ + V^log77|^ dxdt 
Jo Jn 



< E{go,uo,r]o). 

2. Let us now prove that Mg and Af^ are constants of motion. To prove this, we 
make use of Lemma 3.7 and conclude that for any e > there exists i? > such that 



nnB{o,R) 



VSM dx 

By sending 5, /i — > 0, we gather 

Mg{0) ^ rio dx ^ hm / r/sji dx 

. n S,h^O 



Tjs^h dx + 0{e). 



1] dx + 0(e), 



nnB{o,R) 



which implies that Af,, is a constant of motion. 

From the energy estimate (3.33) and Corollary 3.1 it follows that (for a.e t G (0, T)) 



Qh.s'^ dx < C. 



By the requirements of the potential (2.2), we can for any M > determine a radius 
R such that ^{x) > M, for aU \x\ > R. Hence, 



gsji dx = 



gs,h dx 



B{0,R) 



gs,h dx, 



n\B{0,R) 



where 



/ gs.h dx < Af-i / ^gs.h dx < M-^C. 

Jn\B{0,R) Jn\B(0.R.) 



n\B{0,R.} 



Since this holds for any M > 0, wc conclude that Mg{t) is a constant of motion. 



□ 
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4 Large-time Asymptotics 

The analysis in the previous section yields 

Lemma 4.1 Under the hypotheses of Theorem 2.2, we have 

l-T+2 

J |^"lli"(S2) + ll£'l"Hli''2(f2) + ll^?l"llli''3(0) + llk/|u|||LP4(0)rfi== 0, 

lim / |2V,V^+ V^V,$|2 dxdt ==0, (4.1) 



and 

/ / 0''+'' dxdt < C, for allT>l. 

Here, 

o 37 67 

Pi =2, p2 = — — , P3 = — — , P4=4, 

7 + 7 + D 

r 2 1-1 
e = min|-7 - 1, -j. 

Let the sequence {in}^i be such that t„ — !• od and define the sequences: 

u„{t,x) = u{t+tn,x), gn{t,x) ^ g{t+tn,x), ??„ (t , x) = r/(t+t„ , x) , t e (-1, 2), x G f2. 

Observe that, for each fixed n, the triple {gn,'^n-,iln) is a weak solution to particle 
interaction model in the sense of Definition 2.2. 

By virtue of the previous lemma, and bounded energy, we can conclude the exis- 
tence of functions Qs, p, and r]s, such that 

Qn = Q{t + TVi) ^ Qs weakly in {{—1, 2) x SI) 

Pn = p{t + Tn) p weakly in 1, 2) x fl) (4.2) 

Vn = v{t + Tn) Vs weakly in £^((-1, 2), L^{n)). 

Lemma 4.2 Given (4.2) and the hypothesis of Theorem 2.2, rjn — > a.e in (—1, 2) x 
SI, where rjs is given by 

fn exp(— $) dx 

rjs = a, cxp , = ' , 

Jn 10 dx 

which is the unique solution to the problem 



^xTis -?7sV2:$, r]s dx = rio dx. (4.3) 
Jo Jn 



in Jn 

Proof: From the particle equation (1.4), we observe that 

dtVn = -div^(77„u„) + div^; {^/rhi^ + 2V x^/vi) 

in the sense of distributions on (—1, 2) x S7. Hence, 
2 

il<9(?7„||H/-i,i(0) dt <||?7riU„||Li(_i,2;L-i(0)) 

-1 

+ \\nn\\L''-(-i:2;L'^(n))\\^/n^^ + 2V^^||i2(_i_2;L2(o))- 
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Lemma 4.1 tells us that the right-hand side converges to zero as n — ^ oo. Since in 
addition rjn Gb 1, 2; W^{^})), Lemma 3.9 can be applied to conclude that ?7„ — >■ rjs 
a.e in (—1, 2) x fl. In addition, as in (3.39), we find that WxVn ^xV in the sense of 
distributions on (—1,2) x 57. 

Then, as in (3.49), we calculate 



0< / / |V:,V^+ V^V:.$|^ dxdt 

J-lJK 



< lim 



<0, 



for any compact K. The last equality follows directly from (4.1). Hence, we can 
conclude 

Vx^tT = -V^V^E* a.e in (0, T) x n. 



Clearly, this means that also 

V^?/s = -??sVj;$ a.e in (0,T) X 17. 
Next, we let n — > oo in (1.4) and apply Lemma 4.1 and (4.4) to obtain 

r],(t)t dxdt = 0, V0 e Cf=((-1, 2) X n). 



(4.4) 



-1 Jn 



Hence, 77^ is independent of time. Using this, we deduce 



770 dx = lim — 

n— >oo 3 



77,1 dxdt 



lim 



1 Jn 
2 



rin dxdt + - 



-1 JB{0,R) 



rjn dxdt 



rjs dx + lim 



-1 Jn\B(o,R) 
rjn dxdt 



1 JO\B(0,fl.) 



for any ball 5(0, i?) of finite radius R. Now, in view of (2.2), for any M > there 
exists a radius R such that $ > AI. Consequently, 



1 Jn\B{a,R,) 



rjn dxdt < ^ / Vn^ dxdt < M ^D, 



where D is the constant appearing in Corollary 3.1. Letting A/ — > 00, we conclude 

ris dx = / TjQ dx. (4-5) 
1 Jn 

Since (4.3) is linear in 77^ and VxVs^ this concludes the proof. □ 

By taking the limit n — > 00 in the continuity equation (1.2), keeping in mind 
Lemma 4.1, we obtain 

g,(l3t dxdt = 0, V0 G Cr((-1, 2) X n). 



1 Jn 
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Thus, Qs is independent of time. By virtue of Corollary 3.1 and the energy estimate 
(2.6), we have that G L°°{0,T; L^{fl)). Then, we can redo the arguments 

leading to (4.5), with Pg replacing rj, to obtain 

Qs dx ~ / Qo dx. (4-6) 
Jn 

Note that wc do not lose mass in the largc-timc limit (cf [14]). 

Lemma 4.3 Given the convergences (4.2), gi„ — > Qs in L'^ [[—\,2) x fi), where 

+ Q]+V"\ (4.7) 



07 

Cg is uniquely determined by qq dx, and is the unique solution of 

^xP{Qs) = -PQs^x^ in ^, I gs dx = / go dx 

Jn Jn 

Proof: With the already obtained bounds, there is no problem with passing to 
the limit in the momentum equation (1.3) to obtain 

^xP + ^xV = -l^gsS/x"^ - ??sV^$, 

in the sense of distributions on D,. By virtue of the previous lemma, 

V^p= -/3gsS/^^. 

However, then wc arc in the same situation as in [14] and the arguments of [14, 
Proposition 4.1] can be applied to obtain strong convergence p„ gs in L''{{~1, 2) x 
fl). Consequently, p = p{gs), and. in view of (4.6), gs solves 

^xp{gs) = -(igs^x^, I gs dx ^ I go dx 

Jn Jn 

in the sense of distributions. According to [13, Theorem 2.1], gs is then of the form 
(4.7), where Cg is uniquely determined by gs dx. □ 

Lemmas 5.2 and 5.3 concludes the proof of Theorem 2.2. 
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